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Abstract. 

Recently WMAP and BOOMERanG experiments have set stringent constraints 
on the polarization angle of photons propagating in an expanding universe: Aa = 
(— 2.4± 1.9)°. The polarization of the Cosmic Microwave Background radiation (CMB) 
is reviewed in the context of nonlinear electrodynamics (NEED). We compute the 
polarization angle of photons propagating in a cosmological background with planar 
symmetry. For this purpose, we use the Pagels-Tomboulis (PT) Lagrangian density 
describing NEED, which has the form L - {X/K^y-^ X, where X = jF^/^F"^ , and S 
the parameter featuring the non-Maxwellian character of the PT nonlinear description 
of the electromagnetic interaction. After looking at the polarization components in 
the plane orthogonal to the (a;)-direction of propagation of the CMB photons, the 
polarization angle is defined in terms of the eccentricity of the universe, a geometrical 
property whose evolution on cosmic time (from the last scattering surface to the 
present) is constrained by the strength of magnetic fields over extragalactic distances. 



PACS numbers: 98.62.En, 98.80.Es 
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1. Introduction 

Modifications to the standard (Maxwell) electrodynamics were proposed in the 
hterature in order to avoid infinite physical quantities from theoretical descriptions 
of electromagnetic interactions. Born and Infeld [1], for instance, proposed a model in 
which the infinite self energy of point particles (typical of Maxwell's electrodynamics) are 
removed by introducing an upper limit on the electric field strength, and by considering 
the electron as an electric particle with finite radius. Along this line, other models 
of nonlinear electrodynamics (NLED) Lagrangians were proposed by Plebanski, who 
also showed that Born- Infeld model satisfies physically acceptable requirements [2]. 
Consequences of nonlinear electrodynamics have been studied in many contexts, such a, 
for example, cosmological models [3], black holes and wormhole physics |11[5], primordial 
magnetic fields in the Universe [HI [HI |8] , gravitational baryogenesis [8] , and astrophysics 

[niin]. 

In this paper we investigate the CMB polarization of photons described by nonlinear 
electrodynamics. We compute the polarization angle of photons propagating in an 
expanding Universe, by considering in particular cosmological models with planar 
symmetry. The polarization angle does depend on the parameter characterizing the 
nonlinearity of electrodynamics, which will be constrained by making use of the recent 
data from WMAP and BOOMERANG. This kind of investigations has received a lot 
of interest because they represent a probe of models beyond the standard model, which 
may violate the fundamental symmetries such as CPT and Lorentz invariance |l3l[Tl]. In 
what follows we will follow the main lines of the paper on "Cosmological CPT violation, 
baryo/leptogenesis and CMB polarization" by Li-Xia-Li-Zhang [6]. 

2. Minimally coupling gravity to nonlinear electrodynamics 

The action of (nonlinear) electrodynamics coupled minimally to gravity is 

5 = ^ y d'^x^R d^x^L{X, Y) , (1) 

where k = SirG, L is the Lagrangian of nonlinear electrodynamics depending on the 
invariant X = \F^yF^''' = -2(E'^-B^) and Y = lF^^*F^"', where F'''' = V^^A^-V^A^, 
and *F'"' = et'^P'^F^^ is the dual bivector, and e"^^^ = e"'^^^, with e^^^T^ the Levi- 
Civita tensor (£0123 = +!)• 

The equations of motion are [9] 

V.i-LxF^^'' -Ly*Fn = 0, (2) 

where Lx = dL/dX and Ly = dL/dY ^ 

V^F.A + V.Fa^ + VaF^. = . (3) 

After a swift grasp on this set of equations one realizes that is difficult to find 
solutions in closed form of these equations. Therefore to study the effects of nonlinear 
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electrodynamics, we confine ourselves to consider the abelian Pagels-Tomboulis theory 
[16] . proposed as an effective model of low energy QCD. The Lagrangian density of this 
theory involves only the invariant X in the form 

L{X) = - (^) ' ^ = -iX' , (4) 

where 7 (or A) and 5 are free parameters that, with appropriate choice, reproduce 
the well known Lagrangian already studied in the literature. 7 has dimensions 
[energy] 

Following Kunze [9] , the energy momentum tensor corresponding to the Lagrangian 
density L{X) is given by 

T,, = ^ [Lxg'-^F.^Fp, + g^^L] (5) 

and the decomposition of the electromagnetic tensor with respect to a fundamental 
observer with 4- velocity {u^u^ = ~1) 

F^y = 2E[f,u^] - r]f,^^rU'^B^ . (6) 

The electric and magnetic fields are therefore given by Ef^ = F^^u^ and = 

The energy density turns out to be 

1 L 



{26 - 1)E^E'' + B^B" . (7) 



The positivity of p (weak energy condition) imposes, in general, the constraint on 6. 
For the Lagrangian one gets 5 > |. However, this condition can be relaxed because 
we shall consider cosmo logical scenarios in which the electric field is zero, and only the 
magnetic fields survive (this is justified by the fact that during the radiation dominated 
era the plasma effects induce a rapid decay of the electric field, whereas magnetic field 
remains (see the paper by Turner and Widrow in [25])). 

The equation of motion for the Pagels-Tomboulis theory follows from Eq. ([2]) with 
Y = 

V^.F'^'^ = - (5 - 1 ) Ff"" . (8) 
X 

In terms of the potential vector A^^, and imposing the Lorentz gauge V^A'^ = 0, Eq. 
dH]) becomes 

V^V^A^ + R'^^A'' = -{6 - 1)^^ (V'^A'^ - V'A'') , (9) 

where the Ricci tensor R"^ appears because the relation [V^, Vi,]A'^ = —R^^A^. 

To proceed onward, we apply the geometrical optics approximation. This means 
that the scales of variation of the electromagnetic fields are smaller than the cosmological 
scales we consider next. In this approximation, the 4- vector A^[x) can be written as 

m 

A^'ix) = Re [(a^(x) + eb^ix) + . . .)e*^(")/^] (10) 
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with e ^ 1 so that the phase S/e varies faster than the amphtude. By defining the wave 
vector = V^S*, which defines the direction of the photon propagation, one finds that 
the gauge condition imphes /c^a^ = and fc^fo'^ = 0. It turns out to be convenient to 
introduce the normahzed polarization vector so that the vector can be written as 

af'ix) = A{x)e^' , e^.e'' = 1 . (11) 

As a consequence of ( ITTi) . one also finds k^e'^ = 0, i.e. the wave vector is orthogonal to 
the polarization vector. 

By making use of Eq. (|TOl) . one obtains 



^ = ^[l + fi(e)], (12) 

where 



To leading order in e, the term depending on the Ricci tensors can be neglected 
in ([9]). Inserting (fTOl) into ([9]) and collecting all terms proportional to and one 
obtains 



^ : k^k^a" = 2{S - l)A;^A;['^a"l , (13) 

i : 2k^V>'a'' + a'^V^k'' = -2(5 - l)k^V^''a''^ . (14) 

Taking into account the gauge condition fc'^a^ = 0, the first equation implies 

{26-l)P = k^k>' = 0, provided S^^, (15) 

hence photons propagate along null geodesies. Multiplying Eq. f lT^ by a^j, and 
using (fTT]) one obtains 

iv^A;^ = -6k^V^ \nA + {6- 1) V'^V.e^ , 
so that Eq. f|T4|) can be recast in the form 

FV^e" = ^-j^ (16) 

where 



= k^ [V^e" - (e^Vpe'^)^"] 



(17) 
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3. Cosmological setting: Space-time with planar symmetry — )■ universe 
eccentricity — > polarization angle 

3.1. Space-time anisotropy and magnetic energy density evolution 

Let us consider cosmological models with planar symmetry, i.e., having a similar scale 
factor on the first two spatial coordinates. The most general line-element of a geometry 
with plane-symmetry is [20] 

ds^ = dt^ -b'^idx'^ + dy^) - c^dz'^ , (18) 

where b{t) and c(t) are the scale factors, which are normalized in order that 
b(to) = 1 = c(to) at the present time to. As Eq. ( ITSl) shows, the symmetry is on the (xy)- 
plane. The coherent temperature and polarization patterns produced in homogeneous 
but anisotropic cosmological models (Bianchi type with a Friedman-Robertson- Walker 
limit has been studied in [T5]). 

The Christoffel symbols corresponding to the metric ( ITSjl are 

r?, = r°2 = &&, r% = cc, (19) 

pi _ p2 _ ^ p3 _ ^ 
01 ~ 02 ~ ^ ' 03 ~ ^ ■ 

The dot stands for derivative with respect to the cosmic time t. 

To make an estimate on the parameter 6, we have to investigate in more detail the 
geometry with planar symmetry. As pointed out by Campanelli-Cea-Tedesco (CCT) in 
|21] . the most general tensor consistent with the geometry f|T8|) is 

= diag{p, -py, -p\\,-p±) = T^^j^^ + Tj^j^-^^ , 

in which T^j-^^ = diag{p, —p, —p, —p) is the standard isotropic energy-momentum 
tensor describing matter, radiation, or cosmological constant, and = 
diag{p^, —p^, —p^. — p^) represents the anisotropic contribution which induces the 
planar symmetry, and can be given by a uniform magnetic field, a cosmic string, a 
domain wall [22]. In what follows, we shall consider a Universe matter dominated 
{p = 0) with planar symmetry generated by a uniform magnetic field B{t). 

Magnetic fields have been observed in galaxies, galaxy clusters, and extragalactic 
structures [21] , and it is assumed that they may have a primordial origin [23 [H] ■ Due to 
the high conductivity of the primordial plasma, the magnetic field evolves as B{t) ~ 
being frozen into the plasma [23l El] (see below). Denoting with pB the magnetic field 
density, the energy-momentum tensor for a uniform magnetic field can be written as 
^(B) u = pBdiag{l, -1, -1, -1). 

According to (C!) , we find that the energy density of the magnetic field is given by 
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The evolution law of the energy density is given byjlOj 



(21) 



where B is the volume expansion (contraction) scalar, cjab is the shear, and 11"^ the 
anisotropic pressure of the fluid. In a highly conducting medium we still have with good 
approximation S ~ provided that anisotropics can be neglected (this means that 
we neglect radiative effect of the primordial fluid). 

3.2. Space-time eccentricity and polarization angle 

We shall assume that photons propagate along the (positive) x-direction, so that 
= (/c°,fc^,0,0) [7]. Gauge invariance assures that the polarization vector of photons 
has only two independent components, which are orthogonal to the direction of the 
photons motion. Therefore, we are only interested in how the components of the 
polarization vector (2 and 3) change. It then follows that T°" defined in f|T7j) assumes 
the form 



The components of T'^ given by fl22l) vanish in the case of a Friedman- Robert son- 
Walker geometry. 

By defining the affine parameter A which measures the distance along the line- 
element, k^ = dx^/dX, one obtains that e"^ and satisfy the following geodesic equation 




b 



(22) 



(from Eq. (fT6D) 




^- 1 
6 




+ bb{e'f + cc{e'Y 



dX c 



6-1 
6 




+ bb{e'f + cc{e^f 



These equations can be 



further simplified if one observes that fc' 



dt/dX 




(23) 



(24) 



where 



V = k^V^. 



(25) 
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as 



Moreover, the difference of the Hubble expansion rate h/h and c/c can be written 
be 1 de^ 



b c 2(1 - e2) dt 
where we have introduced the eccentricity 



(26) 



The polarization angle a is defined as a = arctan[(c5'^)/ (te^)]. Its time evolution is 
governed by equation 

Pa - ^ (^^ - fj [{bey + ieey] = . (28) 

However, Eqs. fl23|) and fl24|) implies that both be"^ and ce^ evolves as Ai + {S—l)fi(t), 
i = 2,3 , where fi{t) is a function of time and Ai are constant of integration. Therefore, 
to leading order (5 — 1) Eq. (128!) reads 

I)a - ^ ~ c) + ^^^^ - ^^'^ = ° • ^2^^ 

where K = A2 + A^. 

To compute the rotation of the polarization angle, one needs to evaluate a at two 
distinct instants. In the cosmological context that we are considering is assumed that 
the reference time t corresponds to the moment in which photons are emitted from 
the last scattering surface, and the instant to corresponds to the present time. One, 
therefore, gets^ 

Aa = a{t) - a{to) = ——Ke'^iz) , (30) 

4o 

where we have used e(to) = because of the normalization condition &(to) = c(to) = 
1 and log(l — e^) ~ — e^. 

Notice that for 5 = 1 or = there is no rotation of the polarization angle, 
as expected. Moreover, in the case in which photons propagate along the direction z- 
direction, so that k"- = (cuq, 0, 0, k), we find that the NLED have no effects as concerns 
to the rotation of the polarization angle. 

As arises from (!30!) . Aa vanishes in the limit 5 = 1, so that no rotation of the 
polarization angle occurs in the standard electrodynamics, even if the background is 

I Preliminary calculations [40j performed in terms of the electromagnetic field F^^, and of time evolution 
of the Stokes parameters /, Q, U, V (this approach is alternative to one presented in the Sec. II of the 
paper where the analysis is performed in terms of the 4-potential A^) yield again the result (|30p . 
Calculations show that the total flux / is not the same along the three spatial directions, as expected 
owing to the different expansion of the Universe along the x, y and z directions. Moreover the time 
evolution of the Stokes parameters turns out to be a mixture of each others, which reduce to standard 
results as i5 = 1. The polarization angle is defined as 2a = arctan(J7/Q). 
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described by a geometry with planar symmetry. Moreover, even if 5 7^ 1, Aa still 
vanishes for an isotropic and homogeneous cosmology described by the Friedman- 
Robert son- Walker element line (6 = c) ds^ = dt^ — }p'{dx'^ + dy'^ + dz^), because in 
such a case the eccentricity vanishes (this agrees with the fact that for this background 
the components of T*^, Eq. ( |22l) . are zero). 



3.3. Eccentricity evolution on cosmic time 

The time evolution of the eccentricity is determined from the Einstein field equations 
_2_fM + 3//,(ee) + J£j,^2«p«. (31) 

where Hh = b/b. 

It is extremely difficult to exactly solve this equation. We shall therefore assume 
that the e^-terms can be neglected. Since b(t) ~ t^/^ during the matter-dominated era, 
Eq. fl3Tl) implies 



e\z) = 18Fs{z)Q^P , (32) 
where we used 1 + z = b(to)/b{t), e{to) = 0, and 

3 „ 3(1 + ^)^^-3 



(9 - 86) {46 - 3) (9 - 86) {46 - 3) 

is the present energy density ratio 



2{l + z)^. (33) 



,(0) _PB _ B\t,) (B\t,)\-' ( B{to) V fB\t, 



with Per = 3H^{to)/K = 8.1/i2io-47 GeV^ {h = 0.72 is the little-/i constant), and 
i?(to) is the present magnetic field amplitude. 
From Eq. fl30|) then follows 

Aa=^-—^Ke\z,,,). (35) 
46 

where e{zciec)^ the eccentricity ([22]) evaluated at the decoupling z = 1100. 

3.4. Constraints on parameter A from extragalactic B strengths in an ellipsoidal 
Universe 

To make an estimate on the parameter 6, we need the order of amplitude of the present 
magnetic field strength B{to). In this respect, observations indicate that there exist, in 
cluster of galaxies, magnetic fields with field strength (10"'' — 10"^) G on 10 kpc - 1 
Mpc scales, whereas in galaxies of all types and at cosmological distances, the order of 
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magnitude of the magnetic field strength is ~ 10 ^ G on (1-10) kpc scales. The present 
accepted estimations is |2l] [§| 

i?(to)<10-9G. (37) 

Moreover, for an ellipsoidal Universe the eccentricity satisfies the relation < < 
1. The condition > means > 0, with Fs defined in fl33|) . The function F^ given 
by Eq. (133|) is represented in Fig. [H Clearly the allowed region where Fs is positive does 
depend on the redshift z. On the other hand, the condition < 1 poses constraints on 
the magnetic field strength. By requiring < 10^^ (in order that our approximation 
to neglect e^-terms in fl3Tl) holds), from Eqs. fl32|) - fl34|) it follows 

5(to) <9x 10-«G. (38) 

It must also be noted that such magnetic fields does not affect the expansion rate 
of the universe and the CMB fluctuations because the corresponding energy density is 
negligible with respect to the energy density of CMB. 



4. Light propagation in NLED and birefringence 

In this Section we discuss the modification of the light velocity (birefringence effect) for 
the model of nonlinear electrodynamics L{X, Y). We shall follow the paper [31] (see also 
01 [35]), in which is studied the propagation of wave in local nonlinear electrodynamics 
by making use of the Fresnel equation for the wave covectors /c^. The latter are related 

to phase velocity v of the wave propagation by the relation ki = — ki, where ki are the 

V 

components of the unit 3-covector. Thus, in what follows we confine ourselves to the 
phase velocity. It is straightforward to show that for the models under consideration 
(jl]) the group velocity is always greater or equal to the phase velocity [3^ . 
The main result in Ref. [M] corresponds to the optic metric tensors 

= Xg'^- + (y + ,/y - XZ)t^'' , (39) 

(7f = Xg^- + {y- - XZ)t^'' , (40) 

§ The bound (j36p is consistent with the estimation on the present value of the magnetic field strength 
obtained from Big Bang Nucleosynthesis (BBN). As before pointed out, the magnetic fields scales as 
B ^ where the scale factor does depend on the temperature T and on the total number of effectively 
massless degree of freedom g^,s as & (x g^g^^T~^ [28]. The upper bound on the magnetic field at the 
epoch of the BBN is given by [29] B(Tbbn) ^ lO^^G, where according to the standard cosmology 
Tbbn = 10^K~ O.lMeV. Referred to the present value of the magnetic field, the bound on -B(Tbbn) 
becomes [Ml El] 

B{to) = f;^)'s(TBBN) < 6 X 10-^G, (36) 



where Tq = Tito) ~ 2.35 x IQ-'^eV and 5,s(7bbn) g*s{To) ~ 3.91 [2H]- 
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Figure 1. In this plot is represented Fg vs 5 for 5 <\ (upper plot) and 6 >1 (lower 
plot). The condition that the eccentricity is positive follows for Fs > 0. 



which describe the effect of birefringent hght propagation in a generic model for 
nonhnear electrodynamics. The quantities X, y, and Z are related to the derivatives 
of the Lagrangian L{X, Y) with respect to the invariant X and Y, and t^'^ = 

For our model, expressed byEq. (jl]), the quantities X, y, and Z are given by 

X = K^^ = f^X'^'-'\ y = K,K2 = ^{6-l)X^^'-''>-\ Z = 0, 

dL d'^L ^ ^ 

where Ki = 4—— and K2 = 8— —7, while the metrics (13911 and (Mm are 
oX oX^ 

g>r = K,{K,g^^ + 2K2tn , = Klg^\ 

As a consequence, birefringence is present in our model. This means that some 
photons propagate along the standard null rays of spacetime metric g^'^ , whereas other 
photons propagate along rays null with respect to the optical metric Kig^'^ + 2K2t^'^ . 

The velocities of the light wave can be derived by using the light cone equations 
(effective metric) 

Qi'kf.ky = and g^^k^k^ = . 

It is worthwhile to report the general expression for the average value of the velocity 
scalar [3l] 

2. 4 T00(3^ + Ztoo) 2 2 23^2 _ Z + Z(t00)2 + 23^ZtOO 



3A' + 23^too + Z(too)2 3 [X + 2yt'^^ + Z{t 



00^212 
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where T^^ = -t^^ + X = (E^ + B2)/2 {t^^ = -E^), and = S^J^H^"", where 
S = E X B is the energy flux density. The subscript 7 is introduced for distinguishing 
the photon field from the magnetic background. The value of the mean velocity has 
been derived averaging over the directions of propagation and polarization. For our 
model, we get 

+ (41) 

4 r°° 4 

3 4X + 2(5 - l)too ' ~ 3 [4X + 2(5 - l)too]2 

The high accuracy of optical experiments in laboratories requires tiny deviations 
from standard electrodynamics. This condition is satisfied provided |5 — 1| ^ 1. 
Moreover, there are two aspects related to pij) : 

• The average velocity does depend on {only) the parameter 6, so that 7 or A in our 
model can be fixed independently. This task is addressed in the next Section. 

• Because R is positive, one has to demand that 5 — 1 < in order that f ^ < 1. 

The above considerations hold for fiat spacetime, and can be straightforwardly 
generalized to the case of curved space time [M] . 

5. Stokes parameters, rotated CMB spectra and constraints on parameter A 

The propagation of photons can be described in terms of the Stokes parameters /, Q, 
U, and V. The parameters Q and V can be decomposed in gradient-like (G) and a 
curl- like (C) components [32] {G and C are also indicated in literature as E and B), 
and characterize the orthogonal modes of the linear polarization (they depend on the 
axes where the linear polarization are defined, contrarily to the physical observable / 
and V which are independent on the choice of coordinate system). 

The polarization G and C and the temperature (T) are crucial because they 
allow to completely characterize the CMB on the sky. If the Universe is isotropic 
and homogeneous and the electrodynamics is the standard one, then the TC and 
GC cross-correlations power spectrum vanish owing to the absence of the cosmological 
birefringence. In presence of the latter, on the contrary, the polarization vector of each 
photons turns out to be rotated by the angle Act, giving rise to TG and GG correlations. 

Using the expression for the power spectra G{^^ ~ / dk[k'^Axito)AYito)], where 
X,Y = T,G,G and Ax are the polarization perturbations whose time evolution is 
controlled by the Boltzman equation, one can derive the correlation for T, G and G in 
terms of Aa 

G[^^ = Gf^ sin 2Aa , G[^^ = Gf^ cos 2Aa , (42) 

II Notice that in Ref.[3T] the analysis did not include the rotation of the CMB spectra, and in Ref.[32] 
the analysis focused on only the TC and TG modes. Other approximated approaches to discuss the 
rotation angle can be found in Refs. [30ll33] . 
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CiGC 


= 1(67^^-^^ sin4A«, 


(43) 




= Cf^ cos^ 2 Aa + Cp^ sin^ 2 Aa , 


(44) 




= Cf^ cos^ 2Aa + Cp^ sin^ 2 Aa . 


(45) 



The prime indicates the rotated quantities. Notice that the CMB temperature 
power spectrum remains unchanged under the rotation. 

Experimental constraints on Aa have been put from the observation of CMB 
polarization by WMAP and BOOMERanG [HEZI. 

Aa = (-2.4 ± 1.9)° = [-0.00277r, -0.02387r] . (46) 

The combination of Eqs. f l46p and (135 p . and the laboratory constraints |5 — 1| ^ 1 
allow to estimate A. 



5.1. Estimative of K 

To estimate A we shall write 



B = 10"^+^ G 
Fs = 2^3/2 ^ ^ ^100 > 1 . 

The bound (HSl) can be therefore rewritten in the form 



10^^ 10"^ 
^ < \s-l\ < ^ 



where 



A 



A 



KIO 



-6+26 



0.24 X l0-56+2f> 



/GeV^ ^ 



A 



-1 s-i 



The condition |(5 — 1| ^ 1 requires ^ ^ 1. It turns out convenient to set 

^=10", a>0{l). 
From Eqs. (ED]) and fl3T]) it then follows 



A = 10-14+^2 

or equivalently 
Log 



A 



GeV 



J_ -|^ga-2fe+6 

K 



4(5-1) 



GeV. 



-11 



[a - 26 + 6 - LogfsT] . 



(47) 
(48) 

(49) 
(50) 
(51) 
(52) 
(53) 



4(5-1) 

The constant K can now be determined to fix the characteristic scale A. Writing 
A = 10^- GeV, where A^=pi = 19, Aqvt = 16 and Aew = 3 for the Planck, GUT and 
electroweak (EW) scales, respectively, Eq. ( 153|1 yields 



K = io«-2f'+6-C 



^ 4(5 -1)A^ 
C = 4^ — ^ < 1. 



(54) 



14-6/2 

In Fig. |2]is plotted Log(A/GeV) vs K for fixed values of the parameters a, b and 
6 — 1. Similar plots can be derived for GUT and EW scales 
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Figure 2. Kmb K for different values of tire parameter 5 —\, a and h. The parameter 
a is related to the range in which 5—1 varies, i.e. —10^'^^" ^ ^ 1 ^ —10^^^", while 
b parameterizes the magnetic field strength B = 10^^+^ G. The red-shift is 2; = 1100. 
Plot refers to Planck scale A = 10'^^ GeV, with A^^pi = 19. Similar plots can be also 
obtained for GUT (Aqut = 16) and EW (Aew = 3) scales. 



6. Discussion and closing remarks 



In conclusion, in this paper we have calculated, in the framework of the nonlinear 
electrodynamics, the rotation of the polarization angle of photons propagating in a 
Universe with planar symmetry. We have found that the rotation of the polarization 
angle does depend on the parameter 6, which characterizes the degree of nonlinearity of 
the electrodynamics. This parameter can be constrained by making use of recent data 
from WMAP and BOOMERang. Results show that the CMB polarization signature, 
if detected by future CMB observations, would be an important test in favor of models 
going beyond the standard model, including the nonlinear electrodynamics. 

Some comments are in order. In our investigation we have assumed that the planar- 
symmetry is induced by a magnetic field. This is not the unique case. In fact, a planar 
geometry can also be induced by topological defects, such as cosmic string (cs) or domain 
wall [dw) |2T]. In such a case, one has [2T] 



dw 



7 



-Q 



(0) 
dw 



and 



where fl 



[l + z] 



(1 + ^) 



+ 4(1 + ^)^/2-7 



2(1 + z 



,3/2 



(55) 



(56) 



(0) 

{dw,cs) 



are the present energy densities, in units of critical density, of the domain 
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wall and cosmic string. At the decoupling, one obtains 



dw 



and 

e^{Zdec) 



10-^, (5T) 



10-^ , . (5^ 

4 X 10-7 ^ 



The analysis leading to determine the bounds on 5 from CMB polarization goes along 
the line above traced. 

Moreover, a complete analysis of the planar-geometry is required to fix the 
parameter 5. From a side, in our calculations in fact we have assumed that the Universe 
is matter dominated. A more precise calculation should require to use (to solve ([SI])) 
the relation 

1 r 1+^ , 

t = — dz , 59 

Ho Jo ,/njTTzfT^h 

where = 0.3, = 0.7 and Ho = 72 km sec"^ Mpc"^ {z = 1100). From the other 
side, a complete study of the Eq. (13T!) is necessary in order to put stringent constraints 
on the parameter S. 

As closing remark, we would like to point out that the approach to analyze 
the CMB polarization in the context of NLED that we have presented above can 
also be applied to discuss the extreme-scale alignments of quasar polarization vectors 
[38] , a cosmic phenomenon that was discovered by Hutsemekers|36] in the late 1990's, 
who presented paramount evidence for very large-scale coherent orientations of quasar 
polarization vectors (see also Hutsemekers and Lamy [37] )0. As far as the authors 
of the present paper are awared of, the issue has remained as an open cosmological 
connundrum, with a few workers in the field having focused their attention on to those 
intriguing observations. Nonetheless, we quote "en passant" that in a recent paper 
[39] Hutsemekers et al. discussed the possibility of such phenomenon to be understood 
by invoking very light pseudoscalar particles mixing with photons. They claimed that 
the observations of a sample of 355 quasars with significant optical polarization present 
strong evidence that quasar polarization vectors are not randomly oriented over the sky, 
as naturally expected. Those authors suggest that the phenomenon can be understood 
in terms of a cosmological-size effect, where the dichroism and birefringence predicted 
by a mixing between photons and very light pseudoscalar particles within a background 
magnetic field can quahtatively reproduce the observations. They also point out at a 
finding indicating that circular polarization measurements could help constrain their 
mechanism. 

% Hutsemekers and Lamy, and collaborators, have presented, in a long series of papers (not all 
cited here) published over the period 1998 to 2008, a tantamount evidence that the alignment of 
quasar polarization vectors is a factual cosmological enigma deserving to be properly addressed in the 
framework of the standard model of cosmology. The papers quoted here arc intended to call to the 
attention of attentive readers the paramount evidence presenting this cosmic phenomenon. 
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Since cosmic magnetic fields have a typical strength of ~ 10~^ — 10~^ G, on average, 
for a characteristic distance scale of 10-30 Mpc, it is our view that such phenomenology 
can be understood in the framework of a nonlinear description of photon propagation 
(NLED) over cosmic background magnetic fields and the use of a planar symmetry for 
the space-time. Specifically, phenomena involving light propagation as dichroism and 
birefringence can be inscribed on to the framework of Heisenberg-Euler NLED, which 
predicts the occurrence of birefringence on cosmological distance scales. We plan to 
present such analysis in a forthcoming communication |40] . 
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